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ABSTRACT 


Wind  power  generating  turbines  operate  under  constant  as  well  as  rapidly  changing 
conditions.  With  fixed  pitch  blades,  many  wind  turbines  are  allowed  to  operate  regardless  of  wind 
conditions  as  long  as  they  are  able  to  produce  more  electricity  than  it  takes  to  get  them  started. 
However,  the  lifecycle  of  the  turbine  blades  is  often  much  shorter  than  expected  because  of  the 
unsteady  aerodynamic  environment  under  which  they  rotate.  Therefore,  the  National  Renewable 
Energy  Laboratory  (NREL)  has  implemented  a  testing  program  to  determine  the  aerodynamic 
conditions,  and  the  frequency  with  which  they  occur,  which  cause  the  largest  amount  of  fatigue 
on  their  variable  pitch,  three  bladed  downwind  horizontal  axis  wind  turbine  (HAWT).  Different 
techniques  will  be  examined  for  analytically  modeling  the  flow  conditions  with  separation  over  a 
rotating  turbine  blade.  Then,  some  different  techniques  for  implementing  a  feedback  control  loop 
will  be  investigated  to  optimize  tbe  movement  of  the  variable  pitch  blades  on  the  NREL  HAWT. 
The  different  methods  analyzed  will  fall  in  the  two-dimensional,  incompressible  area  with  most 
also  being  for  steady  state  conditions.  The  final  objective  is  to  provide  the  reader  with  a 
background  in  dealing  with  the  aerodynamic  conditions  surrounding  a  rotating  wind  turbine  in  an 
unsteady  aerodynamic  environment. 
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LIST  OF  SYMBOLS 


Coefficient  of  Friction 

Coefficient  of  Pressure 

Shape  factor  of  the  boundary  layer,  H  =  5*/0 

Static  Pressure  (measured  on  the  blade  surface) 

Reference  Static  Pressure 

Static  Pressure  measured  at  the  stagnation  point  on  the  blade's  surface 
Total  Pressure  (measured  from  pitot  probe) 

Dynamic  Pressure 
Perfect  gas  constant 
Reynolds  Number 
Freestream  absolute  temperature 
time 

Streamwise  velocity  component 

Velocity  at  the  edge  of  the  boundary  layer,  corresponds  to  that  from  potential  flow 

Velocity  component  in  the  y-direction 
Local  Freestream  velocity 

Coordinate,  streamwise  direction 
Coordinate  perpendicular  to  the  x-direction 
Thickness  of  boundary  layer  or  shear  layer 
Displacement  thickness  of  boundary  layer 
Momentum  thickness  of  boundary  layer 
Dynamic  viscosity 
Kinematic  viscosity 
Air  density 

Shear  stress  at  the  wall 


INTRODUCTION 


Harnessing  the  power  of  the  wind  as  a  clean  renewable  energy  source  has  been  a  goal  of 
researchers  and  industry  for  many  years.  The  origin  of  windmills  to  utilize  the  wind  to  perform 
work  is  a  centuries  old  concept.  However,  it  is  not  until  recently  that  researchers  have  reoriented 
their  focus  from  the  mechanics  of  building  wind  turbines  to  that  of  optimizing  the  wind  energy 
production  while  increasing  the  blade  lifecycles.  Wind  turbines  are  obviously  placed  in  locations 
where  there  is  a  strong  flow  of  air  on  a  regular  basis.  Any  person  who  has  ever  stood  or  driven  in 
the  areas  where  wind  turbines  operate,  knows  how  rapidly  the  wind  can  change  velocities  and 
directions.  This  unsteady  aerodynamic  environment  leads  to  dynamic  loading  on  turbine 
components,  such  as  generators  and  blades,  "far  in  excess  of  their  design  loads"  [1].  The 
unsteady  aerodynamic  effects  include  "dynamic  inflow,  turbulence,  and  dynamic  stall"  [1]  which 
decrease  the  lifetime  of  major  wind  turbine  components  from  their  20-30  year  design  lifetime  [2] 
quoted  in  [1]. 

The  conditions  under  which  a  wind  turbine  operates  create  stalled  sections  of  the  blade 
along  with  portions  that  have  completely  attached  flow.  In  addition  to  the  normal  statically 
stalled  condition,  the  turbine  blades  also  encounter  a  dynamic  stall  situation.  Static  stall  is  defined 
as  the  loss  of  lift  over  a  wing  due  to  flow  separation.  Dynamic  stall  similarly  results  in  a  loss  of 
lift,  however,  it  is  a  cyclic  phenomenon  with  detachment  of  the  flow  followed  by  reattachment. 
Furthermore,  as  described  by  McCroskey  [3]  quoted  in  [4]  "Dynamic  stall  is  the  phenomenon 
under  which  a  vortex  is  seen  to  travel  from  leading  to  trailing  edge  when  flow  separation  occurs." 
The  existence  of  dynamic  stall  on  a  rotating  turbine  blade  is  evidence  of  aerodynamic  loads  which 
are  exceeding  the  design  parameters.  In  the  past  there  has  been  little  concern  about  the  existence 
of  stalled  portions  of  the  turbine  blade  because  there  were  few  methods  to  alleviate  the  problem 
and  their  costs  were  impractically  high  or  the  method  was  simply  not  implementable.  With  the 
decreasing  costs  and  increasing  speeds  of  microprocessors,  it  is  now  possible  that  control  schemes 
can  be  implemented  to  improve  the  performance  of  turbine  blade  stall. 
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Before  attempting  to  improve  wind  turbine  performance,  it  had  to  be  determined  that 
dynamic  stall  occurred  often  enough  to  warrant  the  expense.  Shipley,  et  al,  [1]  examined  the 
dynamic  stall  phenomenon  on  the  National  Renewable  Energy  Laboratory's  (NREL)  "Combined 
Experiment"  on  a  downwind  horizontal  axis  wind  turbine  (HAWT)  and  statistically  determined 
that  it  occurred  in  over  50%  of  their  test  cycles.  It  is  evident  then,  that  dynamic  stall  occurs  with 
a  frequency  that  justifies  the  expense  of  implementing  a  control  scheme  to  increase  the  lifecycle  of 
turbine  blades  in  addition  to  improving  the  efficiencies  of  wind  turbines.  This  paper  examines 
some  starting  points  for  computer  modeling  the  flow  over  a  turbine  blade  and  progressing  that 
model  towards  use  in  a  feedback  control  scheme.  The  intent  is  to  inform  the  reader  of  the  nature 
of  the  aerodynamic  phenomena  which  control  the  flow  over  the  turbine  blade.  In  order  to 
implement  an  adaptive  control  scheme  it  is  necessary  to  start  with  a  simplified  model  of  the  flow 
and  follow  that  by  more  complex  models  which  ultimately  require  large  numbers  of  empirically 
derived  constants.  However,  it  is  possible  that  a  simplified  model  may  suffice  for  a  reliable 
control  algorithm  which  is  relatively  easy  to  implement. 
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BACKGROUND 


The  freestream  wind  velocity  under  which  turbine  blades  rotate  is  in  an  incompressible 
regime  because  it  is  almost  always  under  Mach  0.1.  This  freestream  velocity  translates  into  a 
reasonable  assumption  that  the  flow  is  incompressible  over  the  turbine  blades.  The  existence  of  a 
constant  density  term  simplifies  the  modeling  of  the  flow  and  is  a  valid  assumption  for  most  wind 
turbines.  C.  Hansen  of  the  University  of  Utah  verified  that  the  errors  created  by  making  the 
incompressible  assumption  are  negligible  when  compared  with  calculations  made  for 
compressible,  isothermal  ideal  gas  [5].  However,  allowing  the  incompressibility  assumption  does 
not  eliminate  the  dynamic  and  transient  events  from  the  complete  model.  Furthermore,  the 
Reynolds  number  over  the  surface  of  turbine  blades  is  generally  greater  than  3x  10  ,  which  means 
that  the  flow  can  be  considered  to  be  turbulent.  However,  the  flow  probably  has  some  regions  of 
laminar  flow  towards  the  leading  edge  of  the  blade.  Also,  for  small  freestream  wind  velocities  the 
flow  could  be  highly  laminar.  This  is  especially  true  for  the  NREL  S809  airfoil  which  is  designed 
to  operate  with  a  maximum  amount  of  laminar  flow. 

Regardless  of  the  type  of  flow  over  the  blade,  there  is  still  the  possibility  of  flow 
separation  and  the  resulting  stalled  condition  when  the  airfoil  experiences  increases  in  angle  of 
attack.  The  reason  for  the  separation  is  that  the  fluid  in  the  boundary  layer  has  low  kinetic 
energy.  Therefore,  it  is  more  difficult  to  move  against  a  pressure  gradient.  However,  a  turbulent 
boundary  layer  has  a  greater  energy  level  and  it  is  therefore  more  difficult  for  the  flow  to  separate 
than  in  a  laminar  boundary  layer.  This  is  an  advantage  of  the  flow  conditions  on  a  wind  turbine. 
The  separation  point  on  a  body  ean  be  determined  through  different  metries.  Ultimately,  the  point 
where  the  shear  stress  goes  to  zero  is  the  point  of  separation.  This  location  can  similarly  be 
denoted  by  3  u/3  y  =  0  which  is  the  shear  stress  without  the  viscosity  term.  Furthermore,  static 
stall  occurs  in  the  existence  of  an  adverse  (positive)  pressure  gradient.  Therefore,  the  static 
pressure  distribution  over  the  turbine  blade  is  of  paramount  interest  in  order  to  obtain  the 
necessary  parameters  to  evaluate  the  aerodynamic  conditions.  However,  since  velocity  and 
pressure  are  interrelated,  the  veloeity  distribution  could  be  measured  or  modeled  directly.  For 
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measurement  on  a  wind  turbine,  pressure  data  is  the  only  realistic  method  because  using  multiple 
hot  wire  or  hot  film  anemometers  on  the  blade  surface  for  direct  velocity  measurement  would 
disrupt  the  flow  too  much.  Furthermore,  the  coefficient  of  pressure  is  usually  the  metric  used  to 
determine  when  dynamic  stall  is  probably  occurring.  Therefore,  when  producing  an  analytical 
model  or  implementing  a  control  scheme  for  a  wind  turbine,  it  would  be  normal  to  work  with  a 
pressure  distribution. 

NREL’s  Combined  Experiment  three-bladed,  variable  pitch,  downwind  horizontal  axis 
wind  turbine  is  a  5.05  meter  radius  machine  with  one  fully  instrumented  untwisted  blade  (see 
Appendix  A  for  basic  dimensions  and  layout).  In  addition  to  the  turbine,  there  is  also  an  upwind 
array  for  measuring  the  upstream  wind  velocity  and  direction  along  with  ambient  air  conditions. 
For  a  complete  discussion  of  the  instrumentation  and  test  setup  see  reference  [6].  The  most 
significant  data  collection  is  from  the  chordwise  static  pressure  ports.  Coupling  this  data  with  that 
from  the  total  pressure  pitot  tubes,  the  velocity  can  be  determined  for  the  same  chordwise 
locations  as  the  static  pressure.  Velocity  is  determined  from  the  following  equation: 

-^  =  ^00  =iP  (1) 


However,  according  to  Shipley,  et  al.  [7]  the  measurement  of  the  total  pressure  is  out  of  phase 
with  the  simultaneous  values  for  static  pressure  because  the  total  pressure  probes  are  0.62  meters 
in  front  of  the  turbine  blade  leading  edge.  However,  when  adjusting  for  the  probe  location  the 
total  pressure  measurement  can  be  used  to  calculate  the  local  velocity  values  as  long  as  the  angle 
of  attack  of  the  instrument  is  below  40*^. 

An  alternate  technique  was  developed  by  Shipley,  et  al.  [7]  to  calculate  the  dynamic 
pressure  independent  of  a  total  pressure  measurement.  This  method  is  based  on  the  knowledge 
that  for  incompressible  and  irrotational  flow,  there  exists  a  stagnation  point  (local  velocity  equals 
zero)  where  Cp  =  10  where 


P-P^ 


Pstagnation  P'>= 


(2) 
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The  results  from  this  method  correlate  well  with  the  calculation  based  on  the  total  and  static 
pressure  measurements  and  corrected  for  probe  location.  However,  with  only  one  dynamic 
pressure  value  for  each  spanwise  location,  there  can  be  no  velocity  distributions  in  the  chordwise 
direction.  Therefore,  measuring  the  total  pressure  maybe  the  only  method  which  will  provide  the 
necessary  velocity  profile  for  implementing  some  control  schemes. 

With  a  pressure  coefficient  determined,  the  occurrence  of  dynamic  stall  can  be  evaluated. 
When  there  is  a  rapid  change  in  local  angle  of  attack,  the  turbine  blade  section  can  surpass  the 
static  stalled  condition  and  enter  a  dynamic  stall.  As  Shipley,  et  al.  [1]  writes,  "Turbulence,  shifts 
in  wind  direction  or  magnitude,  wind  shear,  or  upstream  flow  disturbances  can  alter  the  local 
velocity,  and  hence,  create  blade  angles  of  attack  changes  sufficient  to  drive  dynamic  stall."  This 

rapid  change  in  local  angle  of  attack  will  cause  the  development  of  a  peak  pressure  eoefficient 
near  the  leading  edge  of  the  blade  section.  Therefore,  choosing  a  particular  value  of  Cp  which 

represents  the  occurrence  of  dynamic  stall  is  a  straight  forward  approach  towards  recording  its 
existence.  For  the  Combined  Experiment,  the  blade  airfoil  seetion  was  tested  in  the  Colorado 
State  University  wind  tunnel  and  the  maximum  value  of  suction  recorded  was  =  -5.0  [1]. 

However,  this  testing  was  performed  under  steady  aerodynamic  conditions  in  the  wind  tunnel 
which  are  unlike  those  found  in  normal  wind  turbine  operating  environment.  When  the  analysis 
was  performed  on  the  data  from  NREL's  Combined  Experiment,  a  value  of  Cp  =  -10.0  was 

chosen  as  the  minimum  suction  value  for  measuring  the  occurrence  of  dynamic  stall  [1].  Since 
this  number  is  twice  that  of  the  wind  tunnel  data,  there  is  the  possibility  that  more  dynamic  stall 
events  occur  with  less  suction  on  the  leading  edge  of  the  turbine  blade.  Therefore,  this  criterion 
chosen  for  dynamic  stall  determination  will  necessarily  be  refined  for  the  development  of  an 
optimal  control  scheme. 
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TECHNIQUES  FOR  MODELING  THE  FLOW  ON  A  COMPUTER 

One  place  to  begin  the  analysis  of  modeling  flow  over  a  rotating  turbine  blade  is  to  start 
with  the  basic  laminar,  steady,  incompressible,  two-dimensional  case.  This  is  actually  quite  useful 
when  making  a  complete  model  of  the  flow  since  part  of  the  blade  will  probably  start  with  laminar 
flow  and  continue  through  transition  to  fully  developed  turbulent  flow.  Therefore,  a  good 
solution  would  include  laminar  flow  patched  together  with  transition  and  turbulent  models.  Even 
though  it  is  highly  unlikely  that  the  flow  would  separate  in  the  laminar  region,  to  be  complete,  a 
laminar  flow  separation  algorithm  will  be  discussed.  Laminar  flow  separation  is  relatively 
independent  of  Reynolds  number  as  long  as  the  flow  remains  entirely  laminar.  This  means  that  the 
separation  point  becomes  a  fixed  geometric  point  for  a  specific  body  at  a  set  angle  of  attack  [8]. 
Furthermore,  most  laminar  separation  solutions  rely  on  manipulating  the  boundary  layer  equations 
to  get  a  result. 

Stratford  Criterion 

Stratford  [9]  solved  the  boundary  layer  equations  in  two  parts  by  separating  the  inner  and 
outer  layers.  The  inner  layer  is  dominated  by  the  viscous  forces,  but  they  are  balanced  by  the 
pressure  gradient.  "Therefore,  a  transition  across  the  layer  causes  the  pressure  forces  to  be 
balanced  by  viscous  forces  at  the  wall  and  by  inertia  forces  at  the  outer  edge  of  the  boundary 
layer"  [8].  This  allows  the  effects  of  the  pressure  gradient  to  be  considered  separate  from  the 
effects  of  viscosity  in  the  outer  layer.  This  method  relies  on  the  knowledge  of  potential  flow 
distribution  around  the  body.  The  Stratford  criterion  for  separation  as  modified  by  Curie  and 
Skan  [10]  is 


x^C. 


dC 

p 

y  dx  j 


=  0.0104 


(3) 

(4) 


where  is  the  distance  measured  from  the  maximum  potential  velocity  point  to  laminar 
separation  and  is  the  equivalent  length  upstream  of  the  maximum  potential  velocity  given  by 
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Since  this  method  is  based  on  the  inviscid  potential  flow  solution,  it  is  easy  to  implement  when 
limited  data  is  known  about  the  boundary  layer  growth  over  the  body.  Most  methods,  including 
this  one,  do  not  provide  good  results  at  low  angles  of  attack  when  compared  with  experimental 
results,  but  separation  is  more  of  a  problem  when  the  body  is  under  a  high  angle  of  incidence  to 
the  freestream.  The  Stratford  method  provides  a  relatively  easy  solution  to  the  laminar  separation 
problem,  but  it  is  only  for  a  steady  state  case  to  provide  a  baseline  to  which  an  unsteady  solution 
can  be  compared.  However,  when  attempting  to  examine  unsteady  conditions,  the  equations  are 
the  nonlinear  boundary  layer  equations.  Therefore,  the  ultimate  direction  is  to  integrate  the 
boundary  layer  equations  for  incompressible  flow.  This  method  is  examined  in  the  next  section  of 
the  paper. 

Garner  Criterion 

Now,  since  the  flow  is  most  likely  to  be  turbulent  when  it  undergoes  separation,  an 
algorithm  will  be  examined  to  model  this  condition.  The  Garner  criterion  [1 1]  for  turbulent  flow 

separation  describes  the  skin  friction  as  a  power  law  and  the  following  two  parameters: 

and  r  =  — (7) 
14  dx 


where  x  is  the  distance  along  the  wall  and  n  =  6.  Now,  taking  transition  into  effect  Garner 


developed  the  momentum  and  empirical  equations 

dQ  _1 
dx  6 
dH 


' 

(  13V 

C-r 

H  +  — 

1  7  j_ 

=  g5(H-i.4)r_p  _  00135(7/  -  1.4)] 

dx 


(8) 

(9) 


where  ^  =  ("t  „,/p  u]  )Ree'^"  [8].  Solution  of  the  differential  equations  (8)  and  (9)  yield  a  value  of 
H.  When  H  =  2.6  separation  is  predicted  to  occur.  The  relationship  between  T  and  ©(dHfdx)  is 
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empirically  defined  for  different  values  of  H  in  reference  [11].  The  procedures  for  calculation 
using  the  Garner  criterion  (taken  directly  from  reference  [8])  are  as  follows: 

1 .  Determine  the  development  of  the  laminar  boundary  layer  before  transition. 

2.  Decide  at  which  point  on  the  surface  transition  takes  place,  if  at  all. 

3.  Compute  the  value  of  0  immediately  after  transition.  (The  value  of  0  at  the  position  of 
transition  is  determined  from  Step  1,  and  0  is  continuous  through  transition.) 

4a.  If  transition  occurs  at  or  downstream  from  the  position  of  maximum  velocity,  assume  that 
dHjdx  =  0  there,  i.e.,  that  //  =  1.4  -  (r/0.0135)  immediately  after  transition. 

4b.  If  transition  is  upstream  from  the  point  of  maximum  velocity,  assume  that  H  is  constant 

between  the  positions  of  transition  and  the  maximum  velocity.  The  value  of  0  at  the  position 

of  the  maximum  velocity  is  given  by 

M*  0  =  fM*0l  +  0.007623 

^iruix  V.  ^  A  tr  J 

where  the  subscript  tr  denotes  values  at  the  assumed  position  of  transition  and  the  integral  is 
taken  from  the  position  of  transition  to  the  point  of  the  maximum  velocity  and  k  =  ^H+^. 

5.  Finally,  determine  the  point  of  separation  by  integrating  equations  (8)  and  (9). 

The  Garner  criterion  provides  a  rapid  solution  to  the  turbulent  flow  separation  problem  with 
relatively  reliable  prediction  when  compared  with  experimental  data  [8]. 

Total  Variational  Diminishing  Turbulence  Model 

A  second  turbulent  model  has  been  developed  specifically  for  the  NREL  S809  airfoil 
which  is  the  one  used  in  the  Combined  Experiment.  This  research  is  currently  being  performed  by 
Yang,  et  al.  [12]  at  he  Michigan  Technological  University.  They  are  using  a  spatially  second- 
order  symmetric  TVD  (Total  Variational  Diminishing)  scheme  to  model  the  flow  field  for  the 
NREL  S805  and  S809  airfoils.  Their  model  remains  in  the  steady,  two-dimensional,  viscous, 
incompressible,  turbulent  flow  regime  followed  by  the  previous  model.  However,  their  model 
does  not  predict  when  the  airfoil  will  undergo  separation,  but  rather  models  the  pressure 
distribution  over  the  entire  airfoil.  The  resulting  data  compares  well  with  that  from  Delft 
University  of  Technology  wind  tunnel  experimental  data  except  when  the  airfoil  is  at  a  high  angle 
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of  attack  and  a  large  amount  of  flow  separation  occurs  [13].  Therefore,  Yang,  et  al.  are 
continuing  their  research  into  better  turbulence  models  which  will  allow  prediction  of  flow 
characteristics  beyond  statically  stalled  conditions.  Furthermore,  their  method  still  eliminates  the 
unsteady  terms  from  the  governing  equations.  In  the  absence  of  the  unsteady  terms,  there  is  no 
chance  of  modeling  the  dynamic  stall  characteristics  which  are  found  in  the  field  testing  results  of 
the  Combined  Experiment. 

All  of  the  previous  models  dealt  with  two-dimensional  models  of  airfoils.  Knowing  that 
there  are  still  problems  with  this  level  of  research  makes  it  difficult  to  move  towards  fully  three- 
dimensional  flow  models.  However,  some  models  have  been  developed  which  use  two- 
dimensional  airfoil  data  to  predict  turbine  performance.  According  to  Peter  Tu  [14]  from  the 
Wind  Energy  Program  at  NREL,  the  accuracy  of  the  PROP  and  LSWT  codes  are  poor  when 
involving  high  wind  velocities  and  high  angles  of  attack  on  the  NREL  Combined  Experiment. 

Part  of  the  problem  comes  from  the  fact  that  the  necessary  airfoil  data  from  the  wind  tunnel  tests 
is  not  reliable  in  the  post  stall  region  because  different  tests  have  given  results  which  are  too  high 
or  too  low  when  compared  with  the  Combined  Experiment  results.  Therefore,  if  analytical 
methods  can  be  improved  upon,  such  as  the  previously  discussed  method  by  Yang,  et  al.  [12],  the 
reliability  of  the  PROP  AND  LSWT  codes  could  be  increased.  In  addition,  the  analytical  models 
developed  by  Dini,  et  al.  [15,16]  have  made  progress  "in  modeling  airfoil  stall  and  post-stall 
characteristics  using  an  interactive  boundary-layer  and  trailing-edge  vortex  analysis  method"  [14]. 
In  addition  to  improving  the  airfoil  data,  Tu  also  recommends  developing  an  accurate  tip  loss 
model  for  the  wind  turbine  performance  prediction  codes  [14]. 

One  of  the  major  factors  controlling  the  complexity  of  three-dimensional  analysis  over  that 
of  two-dimensional  is  the  problem  associated  with  the  different  flow  conditions  at  the  tip.  This 
problem  is  aceentuated  when  dealing  with  rotating  propellers  because  boundary  layer  separation 
occurs  at  greater  angles  of  attack  than  when  they  are  stationary  as  in  a  wind  tunnel.  This 
phenomenon  is  explained  by  the  additional  fluid  acceleration  in  the  flow  direction  caused  by 
Coriolis  force  which  in  effect  provides  a  favorable  (negative)  pressure  gradient.  The  Centrifugal 
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forces  cause  less  fluid  to  be  transported  towards  the  blade  tip  from  the  root  than  is  lost  outwards. 
Therefore,  the  boundary  layer  is  thinner  towards  the  blade  tip  which  allows  higher  angles  of  attack 
before  separation  [17],  Muesmann  writes  in  his  thesis  that  transition  occurs  at  a  much  tower  Re 
on  rotating  propellers  when  compared  with  a  stationary  blade  [18].  This  decrease  in  transition  Re 
can  largely  be  attributed  to  centrifugal  force.  Therefore,  obtaining  a  good  post  stall  model  does 
not  necessarily  mean  that  the  results  from  a  turbine  prediction  code  would  perform  any  better  than 
a  two-dimensional  model  unless  there  was  a  good  tip  loss  model.  Furthermore,  the  entire  model 
would  have  to  take  into  account  all  of  the  peculiarities  of  a  rotating  turbine  blade  which  are 
different  from  those  of  the  same  blade  stationary  in  a  wind  tunnel  even  accounting  for  three- 
dimensional  effects.  Therefore,  the  two-dimensional  codes  must  be  improved  before  moving  on 
to  accurate  three-dimensional  rotating  turbine  blade  prediction  models. 

The  two  areas  of  concern  for  two-dimensional  computer  flow  models  are  1)  the 
integration  of  better  turbulence  models  and  2)  the  inclusion  of  unsteady  aerodynamic  terms.  Once 
these  problems  are  addressed,  the  ability  of  a  computational  fluid  dynamics  code  to  model 
dynamic  stall  would  be  greatly  enhanced. 
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TECHNIQUES  FOR  IMPLEMENTING  A  CONTROL  SCHEME 

After  seeing  the  challenges  which  face  the  modeling  of  the  flow  over  a  rotating  turbine 
blade  it  is  hard  to  imagine  establishing  a  control  scheme  to  optimize  the  energy  output  and 
longevity  of  the  components  of  a  wind  power  generator.  However,  there  are  some  techniques 
which  allow  starting  at  a  basic  level  to  determine  which  parameters  of  those  that  govern  the  flow, 
significantly  impact  the  optimization  of  the  control  of  the  turbine  blade. 

Preston  Tubes 

One  a  simple  method  determined  by  Preston  [19]  relies  on  the  fact  that  at  the  point  of 
separation  of  the  flow,  the  skin  friction  becomes  zero.  Therefore,  measurement  of  skin  friction  is 
a  way  to  allow  feedback  to  a  control  loop  which  will  maintain  attached  flow  by  modifying  the 
pitch  of  the  blade.  It  must  be  assumed  that  the  boundary  layer  over  the  blade  is  turbulent  and  the 
region  near  the  surface  follows  the  log  “law  of  the  wall.”  Then,  according  to  Preston  [19]  a  pitot 
tube  can  be  placed  in  contact  with  the  surface  and  wholly  immersed  in  the  law  of  the  wall  region. 
The  difference  between  the  measured  pitot  pressure  P  and  the  measured  static  pressure  p  must 
depend  only  on  p,  v,  and  d  (the  outer  diameter  of  the  tube).  Based  on  dimensional  reasoning, 
a  relation  can  be  developed  such  as 

(10) 

pv’  I  pv’  J 

Experiments  have  been  performed  for  a  flat  plate  whieh  determined  the  function  F  and  the 
resulting  relationship  as 

Young  [20]  writes,  "Because  of  the  small  uncertainties  in  the  values  of  the  constants  involved  the 
method  cannot  be  relied  upon  to  be  in  error  by  less  than  5%  (as  is  the  case  with  the  Clauser  plot 
method)  but  it  has  been  found  to  have  wide  applicability  and  its  simplicity  is  a  major  asset." 
Therefore,  measuring  P,  p,  ,  and  T  and  knowing  p,  d,  can  be  determined.  The 

shear  stress  at  the  wall  is  the  direct  component  in  the  skin  friction  coefficient 
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^  w  _  w 


(12) 


This  method  would  require  experimentation  to  improve  the  empirical  function  F.  Furthermore, 
one  or  more  pitot  tubes  would  have  to  be  placed  on  the  blade  with  telemetry  to  the  control 
computer  inside  the  generator  housing.  The  freestream  velocity  could  be  measured  by  a  pitot 
tube  on  the  generator  housing.  A  temperature  measuring  device  would  also  be  located  here. 
Knowing  the  temperature  and  pressure,  the  density  and  kinematic  viscosity  can  be  determined 


RT 


P  = 


CT 


3/2 


T  +  S 


1^ 

V  =  — 


(13) 


where  C  and  S  are  known  constants  for  the  particular  fluid. 

Using  a  generic  setup  such  as  the  previous  one  would  be  relatively  inexpensive  and  with 
the  small  number  of  parts  would  require  a  low  maintenance.  However,  the  existence  of  a  pitot 
tube  in  the  flow  over  the  turbine  blade  would  disrupt  the  smooth  surface  and  possibly  cause 
premature  separation.  Furthermore,  this  method  would  be  restricted  to  steady  conditions  and 
would  therefore  be  unable  to  detect  dynamic  stalling  conditions  which  occur  rapidly.  Overall, 
using  a  Preston  tube  to  quantify  an  oncoming  stalled  condition  would  be  impractical  since  it 
would  require  empirical  testing  to  check  if  it  is  even  a  valid  method  to  use  with  a  rotating  turbine 
blade.  If  the  method  could  be  validated,  then  it  might  be  used  as  a  baseline  for  other  empirical 
tests  of  control  algorithms. 

Prandtl  Boundary  Layer  Equations 

The  biggest  downfall  of  the  Preston  tube  method  is  that  it  does  not  allow  a  real  time 
control  scheme  with  unsteady  conditions  which  is  the  primary  reason  for  developing  a  variable 
pitch  or  otherwise  controllable  wind  turbine  system.  One  way  to  introduce  the  unsteadiness  back 
into  a  method  is  to  start  with  the  boundary  layer  equations  for  a  flat  plate  in  two-dimensional, 
incompressible,  unsteady  flow  that  were  derived  by  Prandtl  in  the  early  part  of  the  20th  century. 
The  dimensional  quantities  are: 

—  +  .^  =  0  (14) 

dx  dy 
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(15) 


du  du  du  Id  p  d^u 

— — Hu— - \-v— — - t; - l-V-  ^ 

dt  dx  dy  pox  ox 

Boundary  Conditions  y  =  0:  w  =  0,  u  =  0;  y  =  °°:  u=U{x) 

U(x)  is  the  velocity  determined  from  potential  flow 

These  equations  are  derived  by  simplifying  the  two-dimensional  Navier-Stokes  equations  by 

estimating  the  order  of  magnitude  of  each  term.  Equation  (14)  is  the  continuity  relationship  and 

remains  unchanged.  Equation  (15)  comes  from  the  momentum  in  the  x  direction.  To  arrive  at 

this  second  equation  for  the  momentum  of  the  boundary  layer,  it  is  assumed  that  the  change  of  v 

and  p  in  the  y  direction  is  equal  to  zero.  This  is  a  good  approximation  for  boundary  layer  flows 

that  remain  attached  and  do  not  have  a  hypersonic  Mach  number  [21].  Furthermore,  equations 

(14)  to  (16)  remain  applicable  to  curved  surfaces  as  well  as  flat  plates  when  the  surfaces  do  not 


have  sharp  edges  [22]. 

In  order  to  satisfy  equation  (16),  it  is  necessary  to  determine  the  potential  flow  which  is 
the  inviscid  portion  of  equation  (15)  with  no  velocity  gradient  in  the  y  direction.  Also,  u=U  which 
is  the  velocity  at  the  edge  of  the  boundary  layer  determined  by  potential  flow.  This  inviscid  flow 
does  not  come  in  contact  with  the  airfoil  because  of  the  existence  of  the  boundary  layer  and  is 


represented  as; 

dU^^dU^  Idp 

dt  d  X  p  0x 


(17) 


In  order  to  represent  a  steady  flow  case,  3  t  terms  are  set  to  zero  in  equations  (15)  and  (17).  A 

steady  state  flow  condition  reduces  (17)  to  Bernoulli's  equation 

p+\^U'^  =const  (18) 


where  all  quantities  are  dimensional. 

Regardless  of  the  steadiness  of  the  flow,  the  pressure  distribution  is  a  known  quantity. 
Therefore,  the  u  and  v  velocity  terms  are  the  only  unknowns  if  the  flow  is  steady.  In  an  unsteady 
flow,  the  time  derivatives  must  come  from  a  known  function.  Therefore,  (14)  and  (15)  form  a 
simultaneous  system  of  two  equations  and  two  unknowns.  As  previously  noted  one  of  the 
assumptions  for  using  the  simplified  two-dimensional  Prandtl  boundary  layer  equations  is  that  the 
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flow  remains  attached.  Therefore,  this  model  could  allow  the  development  of  a  feedback  control 
loop  that  will  prevent  separation  from  occurring,  but  would  be  unable  to  accurately  model  the 
flow  in  a  post  stall  environment. 

The  latter  method  will  be  a  useful  starting  point  for  developing  a  feedback  control  scheme 
to  optimize  the  performance  and  minimize  loading  from  dynamic  stall  by  modifying  the  blade  pitch 
of  the  NREL  Combined  Experiment  HAWT.  Since  pressure  data  is  readily  available  the  pressure 
coefficient  can  be  calculated  in  real  time  as  a  feedback  element  to  increase  or  decrease  the  blade 
pitch  as  necessary.  However,  one  variable  which  has  yet  to  be  addressed  is  the  yaw  angle  that  the 
turbine  has  with  respect  to  the  oncoming  wind  as  it  changes  directions.  This  is  one  parameter 
which  has  not  been  entrained  into  the  control  model,  but  could  have  some  impact  on  the  control 
algorithm.  However,  once  this  initial  control  algorithm  is  in  place,  the  effect  of  the  yaw  angle  may 
become  more  apparent.  Then,  more  detailed  control  algorithms  may  be  developed.  Furthermore, 
a  yaw  control  loop  could  be  introduced  to  drive  this  angle  to  zero  in  an  attempt  to  eliminate  the 
effect  of  the  yaw  angle. 
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SUMMARY 


Different  techniques  were  examined  for  computer  as  well  as  control  modeling  of  the  flow 
around  two-dimensional  rotating  turbine  blades.  The  primary  concern  lies  in  the  area  of  flow 
separation  which  occurs  in  two  different  modes,  static  and  dynamic  stall.  The  static  stall  reduces 
the  performance  of  the  wind  turbine,  but  does  not  cause  structural  problems.  However,  the 
dynamic  stall  is  one  of  the  primary  drivers  of  a  shortened  lifecycle  of  the  turbine  blades.  As  Bass 
writes,  "The  critical  design  factor  is  the  fatigue  property  of  the  blade  and  the  need  to  keep 
vibratory  stresses  within  acceptable  limits"  [23].  The  real  conditions  existing  in  a  field  test  such  as 
NREL's  Combined  Experiment  may  be  too  complex  to  model  on  the  first  or  even  the  fifth  try,  but 
each  time  the  method  can  improve  on  earlier  shortcomings  and  probably  discover  more  areas  for 
improvement.  Therefore,  picking  the  basic,  very  simplified  models  as  starting  points  with  all  of 
their  deficiencies  is  probably  the  wisest  starting  point  for  development  of  computer  and  control 
models  for  the  NREL  Combined  Experiment  downwind  HAWT.  Then,  if  these  basic  models  are 
successful,  more  advanced  computer  models  and  control  techniques  can  be  used  to  handle  the 
more  complex  flow  problems. 
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Appendix  A 


COMBINED  EXPERIMENT  TEST  SETUP 
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